Multistable transport regimes and conformational changes in molecular quantum dots 
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We analyse non-equilibrium transport properties of a single-state molecular quantum dot coupled 
to a local phonon and contacted by two electrodes. We derive the effective non-equilibrium (Keldysh) 
action for the phonon mode and study the structure of the saddle points, which turn out to be 
symmetric with respect to time inversion. Above a critical electron-phonon coupling A c the effective 
potential for the phonon mode develops two minima in the equilibrium and three minima in the case 
of a finite bias voltage. For strongly interacting Luttinger liquid leads A c = 0. Some implications 
for transport experiments on molecular quantum dots are discussed. 

PACS numbers: 73.63.-b, 71.10.Pm, 73.63.Kv 
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Recently the progress in engineering and production of 
nanoscopic electronic devices reached a new summit as 
experimentalists reported successful contacting of single 
molecules thereby opening the way to build molecular 
quantum dots jj], ||. Since bistable FET-based elements 
play the most important role in modern microelectronics 
a fabrication of such devices on a nanometre scale would 
be a natural way of further progress. The molecules in- 
volved in the electron tunnelling process naturally pos- 
sess elastic degrees of freedom which are bound to re- 
spond to the applied bias voltage in some way; this has 
been dubbed - 'conformational change' (see e.g. 
There have been numerous studies of similar systems 
j|, P though none of them has, to our knowledge, concen- 
trated on the behaviour of the internal degrees of freedom 
in the relevant non-perturbative, non-equilibrium regime. 

To understand the basic physics of the conformational 
change, in this Letter we ignore all the mindboggling 
complexity of real molecular dots Q and employ the sim- 
plest possible Holstein-type model with a single phonon 
mode [0 . We study how this phonon responds to the non- 
equilibrium electron environment. If one applied a 'max- 
imum entropy' principle (such as would be invoked to 
explain the time-reversal breaking in driven macroscopic 
systems), one might expect that the molecule would de- 
form is such a way as to facilitate the current flow. By 
developing the full non-equilibrium theory we shall show 
that such scenario, while admissible in principle, is not 
correct for the quantum system in question. Instead 
we find that in this case a generalised non-equilibrium 
free energy is being minimised and multistable transport 
regimes emerge. 

Let us begin with the Hamiltonian of the system, 

H = H [ik] + H m + Hq + (A + XQ)dU + H t . (1) 

The first two terms on the rhs describe non-interacting 
Fermi seas in the left and right electrodes; tpi, Hi = 
and '02, fJ-2 = —V are the corresponding electron anni- 
hilation operators and chemical potentials, V is the bias 
voltage, see Fig. |l|. The third term describes the phonon 
mode , H Q = (MQ 2 + KQ 2 )/2, where K is the elastic 
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FIG. 1: Schematic representation of the system under in- 
vestigation. A single quantum dot level d coupled to a lo- 
cal boson mode Q can be populated via tunnelling from two 
fermionic reservoirs t[>i and t[>2- The latter are kept at a fixed 
chemical potential difference V. The gate voltage Va enables 
controlling of the dot level population. 

constant and M is the oscillator's mass. The quantum 
dot is described by the fourth term and contains only 
one electronic state, which annihilation and creation op- 
erators are denoted by d and eft, respectively. It can be 
occupied only by one electron (the spin effects are unim- 
portant as long as the Hubbard coupling is weak so that 
the system does not enter the Kondo regime, which we 
assume throughout), which is coupled with a coupling 
constant A to the local phonon mode Q and has the off- 
set energy A. The latter can be controlled by changing 
the gate voltage Vg, see Fig. [l| The last term is re- 
sponsible for the transport through the dot and reads 
H t = ■y[d^(ip 1 + tp 2 ) + H.c.], where 7 is the energy inde- 
pendent tunnelling amplitude (we assume that the leads 
are symmetric to simplify formulas). 

In most cases the molecular vibrational modes are 
slower then the electronic degrees of freedom so that it 
is possible to work in the Born-Oppenheimer (or static) 
approximation. Within this approximation the instant 
value of the current is given by: 
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where T — 2irj 2 po is the resonance width (po being the 
local electron density of states). Outside the static ap- 
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proximation, the phonon coordinate Q still possesses its 
own dynamics. Our goal is to derive the effective action 
for Q. The obvious way to proceed is to eliminate the 
leads and dot variables in the functional integral repre- 
sentation. Since we are interested in the non-equilibrium 
situation, we work in the Keldysh representation, where 
all time integrations are performed along a closed path 
C which consists of the time-ordered branch C_ and the 
ani-time ordered branch C+ M. Generalising the proce- 
dure of Rcf. P| , we define the functional, 



Z[Q] 



(T c e~ iX J c 



(3) 



where Tc stands for the contour ordering operation and 
the average is taken over the exact eigenstate (the steady 
state) of Hamiltonian (|l|) with Q = 0. We set A = (it 
will be restored later) and divide the field Q into two 
Keldysh components scaled by A: Q T = XQ{t on C T ). 
The effective action is then S[Q T ] = 5o[Q T ] + i In Z[Q T ], 
So[Q T ) will be specified shortly. One way to calculate 
the effective action is to take the functional derivative of 
Z[Q T ] respect to, e.g., Q_: 
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-i(T c d(t-+0 + )J(t^e- lX fc dtQ W^),(4) 



t- lies on C_ . If both sides of this equation are divided 
by Z in order to cancel disconnected diagrams introduced 
by Q- being different from Q+, then the rhs becomes the 
time-ordered Green's function of the dot level calculated 
in the presence of Q^{t), 



D—(t,t') = -i(Td(t)rf (?)). 



(5) 



We now implement the non-equilibrium generalisation of 
the standard Born-Oppenheimer approximation by cal- 
culating the above function (i.e. the electronic response 
to the phonon) for static time-independent (but different) 
Q^, restoring the time dependence at the end of calcula- 
tion. By means of fermionic functional integration or via 
the method of Ref. flOl we therefore obtain: 
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-V < uo < 



Combining Eqs. (Q) and (||), one arrives at the following 
expression for the Z-functional: 

lnZ/Z = - J dt J Q dQ- J ^-D—{uj)e w0+ 

+F(Q+) , (7) 

where F{Q+) is a yet unknown function depending only 
on Q+ and the explicit time-dependence of Q T has now 
been restored. This procedure results in the loss of non- 
local in time, dissipative terms, see Ref.B, which are 



accessible via perturbative expansions or a generalised 
Wiener-Hopf method. Such terms are only important 
for tunnelling between different molecular configurations 
and will be addressed in the long version p!l| . In the 
above Zq is responsible for the dynamics of the bosonic 
mode without the electron-phonon coupling and is given 
by the functional integral 
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M(Q 2 _ - Q\) - K(Q 2 _ - Ql) . (9) 



The energy integration in Eq. (f^) results in two contri- 
butions: 



In Z/Z 



i I dt J dQ[n d {Q) + if(Q,Q+)] 



+ F(Q + ). (10) 



The first contribution is purely imaginary and is related 
to the occupation probability of the dot level, 

MQ) = jMcotT^Q/T) + cot- 1 ((Q - V)/T)} . (11) 
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The second contribution is real and is given by 
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Performing analogous computation the field Q+ (or on 
symmetry grounds), one identifies the missing function 
as 

rQ+ 

F(Q+) = - dQ[n d (Q) + if(Q,0)]. (13) 
Jo 

Hence the effective action: 

S = S -iJdt{J® dQ[n d (Q)+if(Q,Q+)] 

' dQ[n d {Q)+if{Q,Q)}}- (14) 

The next step is to minimise the action, which yields the 
system of coupled equations (77 = ±): 

KQ^/X 2 + n d {Qr,) - ivfiQr,, Q- v ) = . (15) 

Let us discuss a special class of solutions which satisfy 
Q = Q- = Q + . Function ( [T2] ) vanishes in this case, so all 
these solutions are real and satisfy the following equation 
(we restored the energy offset A of the dot level) 



KQ/X 2 = n d {Q + A) . 



(16) 



3 




Q/r 

FIG. 2: Effective potential energy of the phonon as a func- 
tion of Q/r for three different gate voltages corresponding 
to A/r = 0, 50, 100 (dotted, solid and dashed lines, respec- 
tively). In the inset: the graphical solution of Eq. (lid). The 
straight line represents the lhs with KY / A 2 = 1. The curves 
show the behaviour of the rid (rhs of Eq. (|l^)) for the same 
A/r as in the main graph. 

We set V = first. It is easy to see that for A < 
and A 3> 1/K there is only one solution of this equa- 
tion, see Fig. |[ The former situation corresponds to 
the dot level being placed below the Fermi edges of the 
leads. In this case the equilibrium position of the phonon 
mode lies below the Fermi energy. For large and posi- 
tive gate voltage A ^> X 2 /K the dot level is effectively 
attracted by the Fermi sea and the dot energy level is 
slightly lowered. Unless the slope of the rhs of Eq. Jl6| ) 
at Q — is smaller than —K/X 2 , which corresponds to 
ttK/X 2 > r _1 , there is a third regime when the line 
—KQ/X 2 crosses n,d(Q + A) three times. In this case 
the effective potential for the boson mode acquires two 
dips, see Fig. |[ and the system becomes bistable. A sim- 
ilar bistability in a different (negative-?/) dot model has 
recently been reported in Ref . Jl2] . The non-equilibrium 
case is even richer. The population probability rid now 
possesses two steps as a function of energy. This opens 
the way to tristable regimes as Eq. ( |l6| ) can have five 
solutions, see Fig. |^. As in the equilibrium situation 
multiple solutions are only possible for sufficiently strong 
electron-phonon coupling, irK/X 2 < r -1 . Remarkably, 
the bistability phenomenon should be visible in the non- 
equilibrium transport even without any adjustments of 
the gate voltage. Starting from the equilibrium system 
with A = and increasing the bias voltage, see Fig. ^, 
one inevitably crosses the bistability regime, which re- 
veals itself as a hysteresis in the I — V diagram. 

For Q_ = Q+, the correction to the effective action, 
(|l4|), can be interpreted as an effective non-equilibrium 
free energy of the phonon mode. We looked for but could 
not find time-inversion breaking minima with Q_ ^ Q+. 
Thus we can confirm that for the system under considera- 
tion the 'quantum Onsager reciprocity principle', put for- 
ward in Ref. is correct. We stress that the Keldysh 
action approach does not require any external conjectures 
to hold and allows, at least in principle, for solutions 
with Q- ^ Q+, which would be at odds with Ref. [fijl. 
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FIG. 3: Flot similar to Fig. ^| but for the non-equilibrium 
case with bias voltage V/T — 50. The dashed, solid and 
dotted lines correspond to A/r = 0, 75, 140, respectively. 
The inset shows the graphical solution of Eq. (hq) for the 
non-equilibrium situation. 




FIG. 4: Current- voltage diagram of a dot with zero off- 
set A = as a function of the applied bias voltage V (cir- 
cles). The current I is calculated via Eq. (jij) after solving 
the Eq. @ and is normalised to e 2 j 2 /2n. The dotted lines 
represent the actual behaviour of the current whereas circles 
in between them are the unstable solutions of the Eq. (|l6|). 
For comparison: solid line is the I — V characteristics of a free 
dot (A = 0) with a finite offset A/r = 50, KT/X 2 = 1. 

Also, even around the time-invertion symmetric minima, 
the (slow) phonon dynamics involves fluctuations of two 
independent but coupled fields Q^. (In equilibrium all 
these complications disappear as Q^i are then simply de- 
coupled.) 

With applications to (single-wall) carbon nanotubes 
in mind, we now briefly discuss the case of interacting 
leads (in equilibrium). The adequate description in low- 
energy sector is then given by the Luttinger liquid (LL) 
model 0, [ll| . To describe the leads we employ the open 
boundary bosonization procedure. In this approach the 
electrons can be thought of as chiral particles living on 
the whole real axis. The tunnelling onto and from the 
dot level takes place at x = in both leads. The Hamil- 
tonian has the simplest form in terms of the Bose fields 
01,2 which describe the collective low-density plasmon 
excitations in the corresponding leads, 

Ho = H [iPi] + H [>p 2 } = ^Jdx [(d x ^) 2 + (d x fo) 2 ] . 
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The local (at x = 0) field operators in this representation 
are then given by p5| , 

^i(t) = {2wa )- 1/2 e lMt)/ ^ . (17) 

This equation contains the Luttinger liquid parame- 
ter g, which is related to the interaction strength Uq 
via g = (1 + 4[/ / 7r ) -1 ^ 2 @, an d the lattice constant 
of the underlying lattice model ao- As a next step we 
substitute ([Tt]) into ([!]) and re-write the dot level opera- 
tors in terms of spin-1/2 operators, S y = — i(<r — d)/2, 
S z = dJd — 1/2. After applying a canonical transforma- 
tion, U = exp[iS z (4>i(0) + <j> 2 (0))/2- s /g], to the Hamilto- 
nian of the system, H 1 = WHU, and introducing new 
fields 4>± — (cf>i ± (\>2)j\f2 we arrive at the following 
Hamiltonian: 

H' = H l<f>±] + QS Z - i7v /2/7ra o ^cos[0_(O)/ v /27] 

where we introduced an additional coupling constant v. 
Transformed into this form the Hamiltonian is closely re- 
lated to the Kondo problem |PiJ. The field Q, which 
is assumed to be adiabatically slow in comparison to 
electronic degrees of freedom, plays the role of the mag- 
netic field applied to the localised spin S. The popula- 
tion probability rid is now related to the magnetisation 
rid = (S z ) + 1/2. Correspondingly, the magnetic suscep- 
tibility gives us the derivative of the rid with respect to 
the boson coordinate Q. 

The special point v = is analogous to the Toulouse 
limit in the Kondo problem. Moreover, at g = 1/3 
the Hamiltonian coincides with that of the four-channel 
Kondo problem, see Ref. Pq| . In order to analyse the 
magnetic susceptibility we proceed along the lines of 
Ref. jl5| . The large time asymptotics of the spin-spin cor- 
relation function turns out to be (TS z (t)S z (0)) = Ct- 2g , 
C being a numerical non-universal constant. Therefore 
the temperature (and hence energy) dependence of the 
magnetic susceptibility \ = drid/dQ is (after the appro- 
priate continuation to imaginary time) : 

X = C 2 [ 1/T dT(S z ( T )S z (0)) ~ d T 2 *- 1 , (19) 
Jo 

where T is the temperature and C± is another non- 
universal constant. The last formula has very important 
consequences for the existence of multistable regimes. In 
the case of strongly interacting electrodes, when g < 1/2, 
X is divergent and the slope of rid is infinite. Hence the 
multistable regimes are always possible when the gate 
voltage is adjusted in appropriate way and there is no 
critical electron-phonon coupling. However, x vanishes as 
the temperature or the relevant energy scale approaches 
zero for weak interactions in the leads, for g > 1/2. So, 



for all 1/2 < g < 1 there is a threshold in A above which 
the multistable regimes emerge. 

To conclude, we have investigated transport properties 
of a quantum dot containing one level coupled to a lo- 
cal phonon mode. In the case of non-interacting leads 
we derived the effective action for the phonon, the anal- 
ysis of which reveals the existence of a bistable regime 
in the equilibrium and a tristable regime for finite bias 
voltage. For strongly interacting leads, the conditions on 
the electron-phonon coupling needed to enter multistable 
regimes are relaxed. One obvious avenue for further re- 
search is the phonon dynamics close to the bistability 
where the dissipation rate changes from Ohmic to super- 
Ohmic at finite bias voltage (llj. The bulk interaction 
effects merit further investigation (so, the g =1/2 model 
is exactly solvable by re-fermionization, even at finite 
V .) The most interesting emerging theoretical question 
is about the viability of time- inversion breaking solutions 
in similar systems, this question is currently open. 
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